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Pracovní úkoly 

1. Ověřte, zda jsou pro dané experimentální uspořádání splněny podmínky platnosti Stokesova 
vzorce pro odpor prostředí při pohybu koule, určete Reynoldsovo číslo. 

2. Změřte dynamickou viskozitu olivového a ricinového oleje Stokesovou metodou. 

3. Pro jednu kapalinu proveďte měření s více typy kuliček. Výsledky porovnejte. 

4. Hustotu skleněných kuliček určete pyknometrickou metodou. 

5. Průměr kuliček d změřte dílenským mikroskopem. 

6. Před měřením zvolte pro každou kapalinu odpovídající velikost kuliček. 

7. Tuto volbu proveďte na základě odhadu splnění podmínek platnosti Stokesova vzorce. Pro 
odhad použijte hodnotu viskozity z tabulek. Pro měření jsou připraveny skleněné kuličky o 
průměru 1 až 4 mm. 

[ 1 ] 

Teorie 

Na kuličku padající ve viskózní kapalině působí tři síly. Kromě tíhové a vztlakové je to ještě 
hydrodynamická síla, kterou lze pro malé v vyjádřit Stokesovým vzorcem 

F x = 67 T7]rv (1) 

kde rj je dynamická viskozita prostředí, kterou chceme zjistit, a r je poloměr kuličky. Vztah 
@ P^tí pouze pro rychlosti charakterizované Reynoldsovým číslem Re « 1, které je definováno 
vztahem 


Re 


2 rpv 
V 


( 2 ) 


kde p je viskozita prostředí. [2] 

Před samotným měřením ověříme velikosti Re pro různé kuličky, hodnotu r změříme dílenským 
mikroskopem, hodnotu rj vyčteme z tabulky. Rychlost vypočítáme tak, že na odměrném válci 
vyznačíme vertikální vzdálenost s a změříme čas t, za který jím kulička projde. Potom víme, že 
v = s/t. Je důležité, aby byla rychlost konstantní, proto budeme měřit spíše v dolní části odměrného 
válce, kde se bude kulička pohybovat už ustálenou mezní rychlostí. Pro vlastní experiment zvolíme 
takové konfigurace, aby byla podmínka Re « 1 splněna. Pro výpočet finální hodnoty budeme 
používat vzorec 


2r 2 (pr -p)g 
V 9v (1 — 2, Ar/R) ’ 

kde R je poloměr odměrného válce, px je hustota kuličky a g je místní tíhové zrychlení. [2] 
Protože se jedná o nepřímé měření, je potřeba odvodit vzorec pro celkovou střední chybu. Ten je 
následující: 
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u(rj) = 0,0925926 


(4) 


(: 9 2 r 6 (p-p T ) 2 AR 2 

+ g 2 r 4 R 2 (p — Pt) 2 (i" — 0, 416667R) 2 At; 2 
+ g 2 r 4 R 2 v 2 (r — 0,416667iž) 2 A/9 2 
+ g 2 r 4 R 2 v 2 (r — 0,416667iž) 2 Ap|. 

+ g 2 r 2 R 2 v 2 (r - 0, 833333i?) 2 Ar 2 
+ r 4 R 2 v 2 (p — pr) 2 (r — 0,416667ií) 2 Ag 2 
\| )/ (v 4 (r — 0,416667-R) 4 ) 

Rychlost kuličky vypočítáme stejně jako v předchozím odstavci. Hustotu kuliček změříme po¬ 
mocí pyknometru a laboratorních vah - nejprve na váhy postavíme suchý, čistý pyknometr i se 
zátkou a zapíšeme si jeho váhu jako mj. Poté ho naplníme rozumným množstvím kuliček (tj. aby 
jich nebylo jen tak akorát na pokrytí dna, ale zároveň ne až k zábrusu), jejichž hustotu měříme, 
zazátkujeme, zvážíme a zapíšeme jako ni 2 - Poté do pyknometru až po okraj přilijeme tekutinu, 
jejíž hustotu známe (v našem případě destilovanou vodu, u které změříme teplotu), zazátkujeme, 
pečlivě otřeme, zvážíme a zapíšeme jako m 3 . Nakonec pyknometr vylijeme/vysypeme, naplníme ho 
pouze destilovanou vodou a zvážíme jako m 4 . Hustota měřených kuliček je potom dána vztahem: 


m 2 - mi 

Pt = P -:-, (4) 

rri 4 — 7713 + rri 2 — mi 

kde p je hustota kapaliny. [3] Protože je měření pvknometrem případ nepřímého měření, je 
potřeba odvodit vzorec pro celkovou střední mezní chybu. Ten je následující: 


u(p T ) 


(mi — m2) 2 (Ap 2 (mi — m 2 + ru.3 


— m.4) 2 + p 2 (Am^ + Am 2 )) + (Am 2 + Am^)p 2 (m 3 — nu ) 2 
(mi — m 2 + 7773 — 7B4) 4 


( 6 ) 


Výsledky měření 

Pomůcky 

• Odměrný válec naplněný olivovým olejem 

• Odměrný válec naplněný ricinovým olejem 

• Dostatečné množství kuliček různých poloměrů 

• Pinzeta 

• Teploměr 

• Pásové měřidlo 

• Posuvné měřidlo 
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• Dílenský mikroskop 

• Stopky 

• Pyknometr 

• Laboratorní váhy 

• Petriho miska a sítko 

Podmínky měření 

teplota 22,9 °C 

tlak 971, 6 hPa 

vlhkost vzduchu 18, 5 % 

místo měření Praha, 

Zavedené veličiny 

l délka, (obecná) 

r poloměr kuličky 

R poloměr válce 

h výška,; vertikální dráha kuličky 
t čas 

A t reakční doba člověka 

g tíhové zrychlení 

SG „specific gravity“ - hustota, látky dělená, hustotou destilované vody o stejné teplotě 
p hustota, 

Pt hustota, kuličky 

pv hustota destilované vody 

po hustota oleje 

T teplota 

rj dynamická viskozita 

Re Reynoldsovo číslo 

a, 6, c parametry litovaných funkcí 
u(X) střední směrodatná odchylka veličiny X 

Hodnoty z externích zdrojů 

Z externích zdrojů jsem získal hodnoty veličin: 

g = 9,83 ±0,01 m/s 2 
Ať = 0,627 ±0,0695 s 


[4] (7) 

[5] (8) 
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Vzhledem k tomu, že reakční doba vždy opožďuje reakci a experiment je navržený tak, aby se 
od sebe zpoždění odečetla, zajímá nás pouze směrodatná odchylka reakční doby. V kombinaci s 
přesností stopek 0, 01 s dostáváme 

u(t) = v/2-0,06952 + 0,01 2 s « 0,1 s (9) 

což budeme odteď považovat za směrodatnou odchylku všech časových měření. 

Dále jsem si zjistil SG ricinového oleje /!/ a hustotu vody /2/ při různých teplotách, zčehož 
jsem vypočítal hustotu ricinového oleje /3/ a vytvořil lineární fit závislosti hustoty ricinového oleje 
na teplotě. Poté jsem z tabulky hustot různých olivových olejů při různých teplotách /4/ vytvořil 
lineární fit pro závislost hustoty průměrného olivového oleje na teplotě /5/. 

Nakonec jsem z tabulek zjistil dynamickou viskozitu olivového a ricinového oleje při různých 
teplotách a opět provedl lineární fit, výsledek zapsal do tabulky /10/. 

Naměřené hodnoty 

•Jako první veličiny jsem změřil teploty obou olejů, průměr odměrných válců, vertikální vz¬ 
dálenost mezi značkami na odměrných válcích, mezi kterými budu měřit čas a z lin. fitu vypočítal 
hustoty obou olejů při naměřené teplotě. Tyto hodnoty jsem zapsal do tabulky /6/. 

Poté jsem pomocí dílenského mikroskopu měřil poloměry jednotlivých kuliček, spočítal jsem 
jejich průměr a směrodatnou odchylku pro každý typ kuličky a výsledky zanesl do tabulky /7/. 
Směrodatná odchylka reflektuje jednak to, jak se liší poloměry jednotlivých kuliček, ale také jejich 
křivost/nepravidelnost. Slouží tedy jako orientační ukazatel, jak přesné budou konečné výsledky. 
Preciznější zohlednění nedokonalostí tvaru kuliček by vyžadovalo delší systematické měření, ideálně 
pro každou hozenou kuličku. 

Následně jsem měřil časy pádu všech kuliček pro výpočet Re. Nejprve jsem experiment pro¬ 
váděl s olivovým olejem, kde jsem výsledek pro každou kuličku změřil třikrát, zjistil průměr a 
směrodatnou odchylku a výsledky zanesl do tabulky /8/. Poté jsem měřil ricinový olej, kde jsem 
kvůli větší časové náročnosti měřil časy pouze jednou pro každou z kuliček; naměřené hodnoty jsou 
zaneseny v tabulce /9/. 

Z výsledků v předchozích tabulkách jsem sestavil tabulku /11/, kde jsou všechny veličiny 
potřebné pro výpočet Reynoldsova čísla. Dynamickou viskozitu jsem vypočítal z lineárního fitu 
v tabulce /10/ a naměřené teploty z tabulky /6/. Následně jsem vypočítal samotné Reynoldsovo 
číslo. 

Vzhledem k podmínce Re « 1 jsem pro samotný pokus zvolil kombinace olivový olej - čirá 
kulička, ricinový olej - čirá kulička a ricinový olej - žlutá kulička. Kombinace ricinový olej - 
tyrkysová kulička sice také splňuje podmínku, ale kvůli časové náročnosti jsem se rozhodl, že tyto 
data už nejsou k přesnému určení zapotřebí. 

Následně jsem použil pyknometr ke změření hustoty kuliček. Postupoval jsem podle metodiky 
nastíněné v sekci Teorie, k výpočtu samotné hustoty jsem použil vzorec ([Š]) a k výpočtu směrodatné 
odchylky (|6]), výsledek jsem zapsal do tabulek /12/ a /13/. 

Nakonec jsem měřil časy padajících kuliček, pětkrát pro každou z konfigurací, a výsledky jsem 
zapsal do tabulky /14/ společně se všemi hodnotami relevantními pro výpočet dynamické viskozity. 

r 


b 



Tu jsem vypočítal společně s její směrodatnou odchylkou pomocí vzorců || a 0 a zapsal do 
tabulky /14/. 

Konečné výsledky mého měření jsou pro konfiguraci čiré kuličky v olivovém oleji: 

r) = 130, 2 ± 4, 9 mPa s, 
pro konfiguraci čiré kuličky v ricinovém oleji: 

rj = 1838 ± 23 mPa s, 

a pro konfiguraci žluté kuličky v ricinovém oleji: 

rj = 1239 ± 40 mPa s. 
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| Ricinový olej lil 



Zdroj 

SG • 10 a 

960 

±5 

při 20°C 

[6] 


955 

±20 

při 25°C 

[7] 


950 

±5 

při 40°C 

[6] 


IVoda 121 

Zdroj 

p [k/m 3 ] 

998,203 

±0,0005 

při 20°C 

[9] 

997,044 

±0,0005 

při 25°C 

[9] 

992,25 

±0,05 

při 40°C 

[10] 


Ricinový olej /3/ z tabulek /!/ a 121 


p [kg/m a 



±5 

při 20°C | 


±20 


1 943| 


při 40°C i 


Lineární regrese p(T) = aT + b 


b= 


-0,75 ±0,04 
973 ±1 



Údaje v tabulce: p [kg/m 3 ], Zdroj: [8] 


11 

17 

23 

29 

35 

45 

55 

65 

75 

1GEV 

911,9 

911,5 

910,6 

910,6 

909,4 

909,2 

908,3 

907,8 

907,1 

2GEV 

912,0 

911,4 

911,0 

910,5 

910,0 

908,8 

908,3 

907,6 

907,1 

3SEV 

911,0 

911,2 

911,0 

910,2 

909,8 

908,4 

908,3 

907,6 

906,9 

4SEV 

912,0 

911,4 

910,7 

910,5 

909,6 

909,2 

908,3 

907,2 

906,0 

5SEV 

911,9 

911,5 

911,0 

910,4 

910,0 

909,0 

908,3 

907,6 

907,0 

6PEV 

912,0 

911,5 

910,9 

909,6 

909,9 

909,2 

908,0 

907,6 

906,8 

7PEV 

911,2 

911,5 

910,1 

909,9 

910,0 

909,0 

908,4 

907,3 

905,8 

8PEV 

911,6 

911,5 

911,3 

909,9 

909,9 

902,0 

908,6 

907,5 

906,4 

9IEV 

911,6 

911,4 

911,3 

910,7 

909,9 

908,4 

908,4 

907,5 

907,1 

10AEV 

912,0 

911,2 

911,2 

909,7 

909,6 

908,9 

908,0 

906,7 

907,3 

1PO 

911,8 

911,7 

910,2 

910,7 

909,3 

908,8 

908,6 

908,0 

906,1 

2PO 

911,9 

911,4 

911,2 

910,5 

909,0 

909,1 

908,1 

908,0 

906,4 

3SO 

911,9 

911,3 

911,0 

910,5 

909,0 

909,1 

908,1 

907,8 

906,4 

4SO 

911,6 

911,4 

911,0 

910,2 

909,0 

909,1 

908,0 

907,8 

906,4 

Průměr 

911,8 

911,4 

910,9 

910,3 

909,6 

908,5 

908,3 

907,6 

906,6 

Směr. odch. 

1,1 

1,5 

1,6 

0,9 

2,5 

1,8 

1,1 

1,3 

1,0 


Olivový olej 151 z tabulky /4/ 

a= -0,08 ±0,003 
_ b- 912,6 ±0,1 


| Oleje /6/ - naměřená data plus výsledek fítů /3/ a /5/ 



T ± 0,5 [°C] 

p [kg/m 3 ] 

u(p) [kg/m a ] 

R ± 1 [mm] 

h ± 2 [mm] 

olivový 

23,0 

911 

5 

572 

149 

ricinový 

22,5 

956 

20 

595 

139 


| Poloměry kuliček lil | 

Typ kuličky 

li ± 0,01 [mn 

1 2 ± 0,01 [mr 

r = V 2 O .2 - li) ± 0,014 [mi 

Ar 2 [mm 2 ] 


[mm] 

rel. odch. 

žlutá, 1. měření 

9,73 

6,89 

1,42 

2,016 

žlutá průměr. 

1,418 

93,40% 

žlutá, 2. měření 

11,70 

8,74 

1,48 

2,190 

žlutá odch. prům. 

1,004 


žlutá, 3. měření 

11,68 

8,97 

1,36 

1,836 

žlutá směr. odch. 

1,325 


tyrkys. 1. měření 

10,79 

8,73 


1,061 

tyrkys, průměr. 

1,022 

93,37% 

tyrkys. 2. měření 

10,68 

8,71 

0,98 

0,970 

tyrkys odch. prům. 

0,723 


tyrkys. 3. měření 

11,31 

9,21 


1,103 

tyrkys, směr. odch. 

0,954 


čirá, 1. měření 

11,13 

9,52 

0,81 

0,648 

1 1 1 1 'I "1 

0,807 

93,34% 

čirá, 2. měření 

10,97 

9,36 

0,81 

0,648 

iisBšsasa 

0,570 


čirá, 3. měření 

10,94 

9,32 

0,81 

0,656 

S^ESíEŠSSKSS 

0,753 



| Časy - olivový olej 181 | 

Olivový olej 

t± 0,1 [s] 


[mm] 

rel. odch. 

žlutá, 1. měření 

2,02 

žlutá průměr. 

2,157 

5 , 20 % 

žlutá, 2. měření 

2,27 

žlutá odch. prům. 

0,073 


žlutá, 3. měření 

2,18 

žlutá směr. odch. 

0,112 


tyrkys. 1. měření 

3,36 

tyrkys, průměr. 



tyrkys. 2. měření 

3,35 

tyrkys odch. prům. 



tyrkys. 3. měření 

3,36 

tyrkys, směr. odch. 



čirá, 1. měření 

11,13 

čirá průměr. 


0 , 88 % 

čirá, 2. měření 

10,97 

čirá odch. prům. 



čirá, 3. měření 

10,94 

čirá směr. odch. 

0,097 
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| Časy - ricinový olej 79/ | 

Typ kuličky 

t ± 0,1 [s] 

rel. odch. 

žlutá 

25,38 

0,39% 

tyrkys. 

39,56 

0,25% 

čirá 

56,31 

0,18% 


lViskozita olejů z tabulek 710/ Zdroj: [ll] | 

1 



T [°C] 


Dyn. visk. 

25 

30 

50 

75 

100 

olivový i} = 

67 

54 

25,8 

9,4 

7 

ricinový r| = 

700 

451 

125 

42 

16,9 

1 



f(x) = exp(ax + b) + c | 


a 

b 

c 

Aa 

Ab 

Ac 

olivový 

-0,043856 

5,23935 

3,84088 

0,00415 

0,09153 

1.815 

ricinový 

-0,083546 

8,58785 

27,7702 

0,008458 

0,2118 

13,69 

1 



n 

An 


olivový při naměřené teplotě 

72,60 

6,97 

ricinový při naměřené teplotě 

847 

160 


Reynoldsovo číslo 711/ 

Re = 2rph/tri 

Kulič., olei 

t±0,l [s] 

h ± 2 [mm] 

r [mm] 

u(r) 

p [kg/m“] 

u(p) 

p [mPa s] 

u(n) 

Re 

u(Re) 


2,16 

149 

1,42 

0,05 

911 

5 

72,60 

6,97 

2,46 

0,28 



149 

1,02 

0,03 

911 

5 

72,60 

6,97 

1,14 

0,12 


5,86 

149 

0,81 

0,00 

911 

5 

72,60 

6,97 

0,51 

0,05 

žlut., ricin. 

25,38 

139 

1,42 

0,05 

956 

20 

847 


0,0175 

0,0034 

E3SSH1S 

39,56 

139 

1,02 

0,03 

956 

20 

847 


0,0081 

0,0016 

1133391 


139 

0,81 

0,00 


20 

847 

160 

0,0045 

0,0009 


|Pyknometr — čiré kuličky 712/ | 


m [g] 

u(m) [g] 

prázdný 

9,1997 

0,0005 

s kuličkami* 

19,8989 

0,0005 

s kuličkami a dest. H 2 0 

26,8282 

0,0005 

s dest. H?0 

19,3451 

0,0005 

BBBBzsaaaiígssii 

997,7 

0,2 

Ihust. kulič.: p(T) [kg/m 3 ] 

3319 

31 


Pyknometr - žluté kuličky 713/ i 


m [g] 

u(m) [g] 

prázdný 

9,3278 

0,0005 

s kuličkami 

18,6617 

0,0005 

s kuličkami a dest. H 2 0 

25,0070 

0,0005 

s dest. H 2 0 

19,3701 

■H 



0,2 

Ihust. kulič.: p(T) [kg/m 3 ] 

mmn 

24 


Naměřená viskozita ole 

ů 714/ 

Konfigurace 

r [mm] 

t±0,l [s] 

P(T) [kg/m 3 ] 

p(O) [kg/m 3 ] 

h ± 2 [mm] 

R ± 0,5 [mm] 

H [mPa s] 

u(n) 

olivový, čirá 

±0,005 


±31 

±5 





0,80 

5,51 

3319 


149 

28,5 

131,8 

3,5 

0,80 

5,59 

3319 

911 

149 

28,5 


3,5 

0,80 

5,44 

3319 

911 

149 

28,5 

131,8 

3,5 

0,75 

5,91 

3319 

HD 

149 

28,5 


3,6 

0,79 

5,58 

3319 

mamm 

149 

28,5 

130,0 

3,5 

ricinový, čirá 

±0,005 


±31 

±20 





0,80 

72,12 

3319 

956 

139 


1808 

39 

0,80 

72,07 

3319 

956 

139 

30 

1830 

39 

0,81 

72,25 

3319 

956 

139 

30 

1882 

40 

0,82 

69,74 

3319 

956 

139 

30 

1840 

40 

0,82 

69,31 

3319 

956 

139 

30 

1829 

39 

ricinový, žlutá 

±0,05 


±24 

±20 





1,47 


2519 

956 

139 


1232 

33 

1,49 

21,87 

2519 

956 

139 

60 

1259 

34 

1,44 

23,07 

2519 

956 

139 

60 

1238 

33 

1,53 

«a 

2519 

956 

139 

WZ 

1274 

34 

1,50 

mmsm 

2519 

956 

139 


1194 

32 
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Diskuze výsledků 

Pokud porovnáme výsledek naměřený při konfiguraci žlutých kuliček v ricinovém oleji, 1239 ± 
40 mPa s, s dynamickou viskozitou ricinového oleje z tabulek 847 ± 160 mPa s /10/, zjistíme, že 
náš výsledek má relativní chybu 1,46. Taková chyba je poměrně vysoká, ovšem s přihlédnutím k 
nejistotám obou hodnot a k variabilitě olejů dostupných na trhu přijatelná. 

Zbylé dvě konfigurace oproti tomu mají velmi vysokou relativní chybu (konfigurace čirá—ricinový 
2,17, konfigurace čirá—olivový 1,81), která naznačuje, že při experimentu došlo k systematické 
chybě. Důležitým vodítkem zde je, že hodnoty pro žluté kuličky zhruba odpovídají realitě, kdežto 
hodnoty pro čiré kuličky jsou hrubě zkreslené. Jedinou měřenou hodnotou, která byla měřena pou¬ 
ze jednou (tj. je náchylná na sytematickou chybu) a přímo souvisí s typem kuličky je pr- Pokud 
tedy zavrhneme naměřenou hustotu čirých kuliček jako nesprávnou a budeme předpokládat, že 
byly oba druhy kuliček vyrobené ze stejného materiálu a mají tudíž i stejnou hustotu, vyjdou nám 
poněkud jiné hodnoty - pro olivový olej to je: 

rj = 86, 9 ± 3,4 mPa s, 

a pro ricinový: 

r) = 1216 ± 18 mPa s. 

Nová hodnota dynamické viskozity pro olivový olej 86, 9 ± 3,4 mPa s velice přesně souhlasí s 
viskozitou získanou z tabulek 72,6 ±7,0 mPa s, relativní chyba výsledku je 1,20. Nová hodnota 
pro ricinový olej navíc souhlasí s hodnotou naměřenou u žlutých kuliček - relativní chyba nové 
hodnoty vztažené k hodnotě žlutých kuliček je 0, 98. 

Tento výsledek je dostatečně přesvědčivý na to, abych považoval systematickou chybu v měření 
čirých kuliček za téměř jistou. Nyní se ještě můžeme zamyslet nad původem této chyby. Naměřená 
hustota kuliček 3319 kg m -3 je o 800 kg m -3 vyšší, než jakou předpokládáme, a závisí na pěti 
parametrech: hustotě vody, hmotnosti prázdného pyknometru, hmotnosti pyknometru s kuličkami, 
hmotnosti pyknometru s kuličkami a vodou, a hmotnosti pyknometru naplněného jen vodou. 

Hustota vody s nejvyšší pravděpodobností příčinou chyby není, protože je vyčtená z tabulek 
a stejná jako u měření hustoty žlutých kuliček. Hmotnost prázndého pyknometru bychom museli 
snížit o 2 g, abychom se dostali na požadovanou hodnotu - není ovšem příliš pravděpodobné, že 
bych na vahách při měření měl pyknometr a navíc ještě nějaký dvougramový předmět. 

Stejně tak by musela být o gram nižší hmotnost pyknometru s kuličkami a vodou. Jediný 
faktor, který by mohl zvýšit naměřenou hmotnost oproti ideálnímu případu, je možnost, že jsem 
před vážením pyknometr dostatečně neutřel. Ale 1 g vody na povrchu pyknometru je přeci jen 
poměrně vysoké číslo. 

Hustota pyknometru naplněného vodou by musela být ve skutečnosti o gram vyšší než naměřená 
hodnota, abychom dosáhli požadovaného výsledku, což nepravděpodobné. 

Nakonec hustota pyknometru s kuličkami by musela ve skutečnosti být o 2 g vyšší než naměřená 
hodnota. To by mohlo být způsobené tím, že jsem při měření pyknometr zapomněl zazátkovat. 
Hmotnost 2 g by rozměrům a materiálu zátky odpovídala. Navíc při tomto měření není v pykno¬ 
metru žádná kapalina, takže je snadné na zátku zapomenout. Proto je tato hypotéza nejpravdě¬ 
podobnějším vysvětlením systematické chyby výsledku. 
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Závěr 


První výsledky hrubě neodpovídaly realitě. To vedlo k objevení systematické chyby v měření 
hustoty čirých kuliček. Po opravě této systematické chyby již výsledky vycházejí správně s vysokou 
přesností. Abychom si však mohli výsledkem být jistí, bylo by třeba hustotu čirých kuliček znovu 
přeměřit. 
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